Theory for transport through a single magnetic molecule: Endohedral N@C 6 o 
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We consider transport through a single NQC60 molecule, weakly coupled to metallic leads. Employing a 
density-matrix formalism we derive rate equations for the occupation probabilities of many-particle states of 
the molecule. We calculate the current-voltage characteristics and the differential conductance for N@C6o in 
a break junction. Our results reveal Coulomb-blockade behavior as well as a fine structure of the Coulomb- 
blockade peaks due to the exchange coupling of the Ceo spin to the spin of the encapsulated nitrogen atom. 

PACS numbers: 85.65.+h, 75.50.Xx 
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I. INTRODUCTION 

The rapid progress of miniaturization of electronic de- 
vices has led to chip features smaller than 100 rim, for 
which standard semiconductor technology reaches its limit. 
One proposed solution is a transistor consisting of a single 
molecule. In recent years transport through single molecules 
has been studied quite extensivelyi^iiii^ — for example, in 
break junctions. 7 ' 8 ' 9 Inelastic transport occurs due to the inter- 
action of electrons with internal vibrational or magnetic de- 
grees of freedom of the molecules. Transport through mag- 
netic molecules 101112 is particularly interesting also from the 
point of view of spintronics — i.e., the idea of exploiting the 
electron spin in electronic devices. While most molecules 
are normally nonmagnetic, there are exceptions such as en- 
dohedral NQCgo — i.e., a nitrogen atom encapsulated in a Ceo 
cageJ^ It is known that the encapsulated atom retains its p 
electrons, 13 leading to a localized spin Sn = 3/2. There 
are fascinating ideas of employing this spin in a quantum 
computerJl 

In this paper we propose to measure the current through 
a single NQCgo molecule in a break junction and we calcu- 
late the current-voltage (I-V) characteristics and the differen- 
tial conductance dl/dV. Since transport through a single Cgo 
molecule has been demonstrated 8 and the synthesis of endohe- 
dral fullerenes is also feasible, 13 such an experiment is possi- 
ble with present-day apparatus. Besides the typical Coulomb- 
blockade behavior we predict a characteristic fine structure of 
the Coulomb-blockade peaks in dl/dV due to the exchange 
coupling of the Cgo spin to the spin 3/2 of the encapsulated 
nitrogen atom. It should be mentioned that the discussion 
of transport through P@Cgo proceeds quite analogously and 
yields qualitatively identical results. 



H. THEORY 

In our model the N@C6o molecule is treated as a quantum 
dot and the leads, labeled as L (left) and R (right), as macro- 
scopic charge reservoirs. Relaxation in the leads is assumed 
to be sufficiently fast so that the electron distributions in the 
leads can be described by Fermi functions. As Ceo generally 
prefers to be singly or doubly negatively chargedjA^ii&iiA&ISi 



we assume that electronic transport through the molecule in- 
volves only the threefold-degenerate LUMO (lowest unoccu- 
pied molecular orbital). 20 Since we concentrate on the fine 
structure of the differential conductance close to degeneracy 
points, we assume that the fivefold-degenerate HOMO (high- 
est occupied molecular orbital), which lies about 7.5 eV below 
the LUMO, 18 remains fully occupied, whereas the threefold- 
degenerate LUMO+1, about 1.7 eV above the LUMO^ re- 
mains empty. 21 Charge transfer from the nitrogen atom to 
the fullerene cage is assumed to be negligible. When the 
LUMO is partially occupied, the net spin of the electrons 
in the LUMO, Sc 6Q , couples to the spin 3/2 of the nitrogen 
atom, Sn, and the total spin of the molecule is S = Sc 60 +Sn- 
The exchange interaction can be written in the simple form 
— JSc 6Q • Sn due to Hund's first rule: Note that the three 
LUMO's and the three nitrogen p orbitals both have odd par- 
ity with a single nodal plane each, which can be chosen as 
the xy, yz, or zx plane. Consequently, there is only a sig- 
nificant exchange interaction between any LUMO and the p 
orbital of the same symmetry. The exchange interaction can 
thus be written as a sum of terms for the three pairs of orbitals. 
However, due to the strong Hund's-rule coupling, the p orbital 
spins combine to Sn = 3/2; i.e., they are all parallel. Then 
the spin of each p orbital is Sn /3, as can be proved by project- 
ing the spins onto the Sn = 3/2 subspace. Thus the exchange 
terms can be combined to the simple scalar product. The full 
Hamiltonian of the system then reads H = H& + iZi ea ds + H t , 
where 

H d = (e - eU g ) n d + ^n A {n A - 1) - J Sc 60 ■ S N (1) 
represents the molecular quantum dot, 
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(2) 
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represents the leads, and 

H X = Y Y( taa aka C no- + C c L a ak ( x) (3) 
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describes the tunneling between the dot and the leads, which 
is assumed to be weak compared to typical excitation ener- 
gies of the molecule. Here, the operator a nrT creates an elec- 
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tron with spin a in the molecular orbital n, while a^ kcr cre- 
ates an electron in lead a with spin a, momentum k, and en- 
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and 



ergy e Q k relative to the Fermi energy. rid — 
s c 60 = Y,naa' 4<r (°W / 2 ) c na> are the number and spin op- 
erators of electrons on the dot, respectively. Electron-electron 
interaction is taken into account by the local Coulomb repul- 
sion U and the exchange interaction with the nitrogen spin by 
the exchange coupling J. The values of e, U, and J are not 
well known at present. Ab initio calculationsi^iii^ indicate 
that Cgo is the ground state, whereas Cgo 2 has a slightly 
higher energy. This is in agreement with the experimental 
observation of a very-long-lived metastable Cgo 2- (Ref. fT^) . 
However, there are recent contradicting ab initio results pre- 
dicting Ceo 2 to be slightly bound relative to Ceo - (Ref.ll^). 
For our numerical calculations we use e = — 2.75eV and 
U = 2.84 eV in accordance with Ref.[nl J appears to be fer- 
romagnetic. We take J ~ 1 meV from ab initio calculations 
of Udvardi. 22 This relatively strong exchange coupling is con- 
sistent with the absence of electron-paramagnetic-resonance 
(EPR) signals for N@Ceo anions with charges —1 through 
—5, while the signal is present for the neutral molecule and the 
hexa-anion. 23 Note that the exchange coupling is significantly 
smaller than the energy of relevant vibrational modes. The 
oscillations of the molecule as a whole have an energy of the 
order of 5 mcV. 8 The oscillations of the nitrogen atom within 
the Ceo have an energy of 13 meV, 24 whereas the modes of 
the Cgo cage lie at much higher energies. 

We next derive rate equations for this model starting 
from the equation of motion for the density matrix p£&& 
dpi(t)/dt = —i [Ha, pi(t)]. Here, the index / denotes the 
interaction representation with respect to H t . Integration and 
iteration give 
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dpijt) 
dt 



= -i[J3tj(t),Pi(0)] 



(4) 



dt' [H t i(t),[H tI (t'), Pl {t') 



Assuming that the leads are weakly affected by the quantum 
dot and neglecting correlations between the two, pi(t) can be 
replaced by the direct product of the reduced density matrix 
of the dot, pa(t) = Tri ea( j s pi(t), and the density matrix pi ea ds 
of the leads, piit) ~ pa it) <8> Aeads- We then obtain 



- / dt 1 Tr leads [H u (t), [H tI (t'), p dI (t') ® p lead s] 
Jo 



dpdijt) 
dt 

(5) 

Returning to the Schrodinger representation and using the 
Markov approximation^ 2 ^ pn(t') « puit), we find 



dpd(t) 



dt 



-i [H d , p^ - Tr^ds / dt' 
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(6) 



as the equation of motion for pd- This expression shows that 
the tunneling Hamiltonian H t is treated in second-order per- 
turbation theory. Taking the trace over the degrees of freedom 
of the leads produces Fermi functions according to 
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Saa'Skk'Saa' f{tak — Pa), (7) 



where p a denotes the chemical potential of lead a due to the 
applied source-drain voltage V. Expanding the nested com- 
mutators in Eq. l|6} and inserting Eq. l|3} gives eight terms: 



dpdjt) 
dt 



/ dt' V \t a 



x {/(e« k - PaV^'c^e^'cl^ 'p A (t) 

+ [1 - /(e«k - Pa)] e-^'c^e-^'^e^'piit) 

- [1 - /(e Q k - Pa)] e 1 ^ 1 ' c naPA {t)e- lH ^' cl^ 1 ' 

- /(e«k - p a )e~ l ^ kt 'ci aPi {t)e- iHit ' c n .„e iH * f 



[1 - /(e Q k - p a )] e 



'** c na e iH ^ pi(t)cl c 



- /(e ak - PaY^ 1 ' e~ lHit ' cX^ 1 ' P d(t)c nla 
+ /(e«k - Pa)e- l ^ t 'p A (t)e-^ t 'c na e^ t 'cl iT 

+ [1 - /(eak - Pa)] e^ 1 ' p A (t)e- lH ^' cl^'cn^). 

(8) 

The probability of the dot being in the many -particle state \n) 
is P n = (n\p A (t)\n). Introducing the overlap matrix elements 

Cmn = H Ei c i*\ n ) and C$ n = (m| cjjra) and iden- 
tifying the integrals in Eq. (JHJ as <5 functions we can write 
Eq. (|8} as a set of rate equations 



dP n 
~dt 



P m E 

1 1 L r . 
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(9) 



with transition rates 



27r|t Q | 2 D Q /(e^ 
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Here, is the energy of the many-particle state |n) of the dot 
and D a denotes the density of states per spin species in lead 
a, which we take to be constant and equal for both leads. The 
matrix elements C^ nn (C^ n ) can only be finite if the electron 
number of state is larger (smaller) by 1 than the electron 
number of state \m). We are interested in the stationary state, 
which corresponds to dP n /dt — for all states \n). 

In deriving Eq. we have assumed that the density ma- 
trix pi is completely diagonal. This assumption requires some 
thought since many of the eigenstates of our molecular quan- 
tum dot are degenerate so that one might expect finite off- 
diagonal components even in the stationary state. However, 
this is not the case: Let U be a unitary matrix that leaves 
the dot Hamiltonian H A invariant. With any stationary den- 
sity matrix p A , UpdU^ is another solution. Now suppose that 
there exists a stationary solution p A that is not diagonal within 
a block of degenerate states. Then one can choose U so 
that UpdW is diagonal since the nonzero off-diagonal com- 
ponents have been assumed to connect degenerate states (we 
exclude the case of accidental degeneracy). But then U pdW 
has unequal diagonal components — i.e., probabilities P n — 
for symmetry-related states. This is clearly unphysical. On 
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the other hand, if is already diagonal with degenerate dot 
states having equal diagonal components, any allowed trans- 
formation U leaves invariant. 

The current operator for lead a — L, R reads 26 

I a = i [H,N a ] = -iy^X t <x C n<T a >otea ~ C a ikcr C no-): C 11 ) 
nk(T 

where N a is the number operator for electrons in lead a. Trac- 
ing out the leads we arrive at an expression for the expectation 
value of the current: 



(I a ) = 2ixD a \t a 



i 2 E 

ml a 



[1-/(4 



4 \cu 2 )p m - 

t R and fj, L 



(12) 



We here consider the symmetric case 
-MR = V/2. 

As there are (1) possible ways of filling the threefold- 
degenerate Ceo LUMO with i electrons according to the Pauli 
principle and as the nitrogen atom possesses a spin 3/2, solv- 
ing the rate equations and calculating the current reduces to 
an eigenvalue problem of dimension 4 x Y^=o (i) = 256. 
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FIG. 1: Current I and average occupation (n) of the Ceo LUMO as a 
function of the source-drain voltage V = fiL — Mr for e = —2.75 eV, 
U = 2.84 eV, J = 1 meV, V g = V, and T = 0.01 K. The inset 
shows the fine structure of one particular Coulomb-blockade step. 



III. RESULTS AND DISCUSSION 

The I-V characteristics plotted in Fig.^show a conductance 
gap for \ V\ < 0.18 V due to Coulomb blockade. Our numer- 
ical results show that the current / is symmetric with respect 
to the applied source-drain voltage V in accordance with the 
high symmetry of the fullerene molecule. Each step in the 
main I-V curve corresponds to the opening of additional cur- 
rent channels. Simultaneously, the average occupation (n) of 
the dot changes. For the parameters chosen above, the Cgo - 
state is the ground stated At the first step, the potential drop 
becomes large enough to allow transitions between the nearly 
degenerate charge states —1 and —2, as the chemical poten- 
tial ml = V/2 reaches the value assumed for the ionization 
energy £(C 60 2 ~) - #(C 60 ~) = e + U = 0.09 eV. At the 
second step, transitions between the charge states —1 and 
become possible, etc. 

Our results for the occupation probabilities reveal that de- 
tailed balance is satisfied for the broad plateaus in Fig.^ — i-e-, 
P n R n ^ m — P m R m ^ n . As a consequence, the dot occupa- 
tion probabilities P n for all occupied states are identical in 
the limit T — > 0, as the transition rates R n ^ m are symmetric 
for each pair n, m of occupied states. This also accounts for 
the fact that the average occupation (n) is exactly unity for 
V = V, increases to (24 x 1 + 60 x 2)/(24 + 60) = 12/7 w 
1.71 at the first step, when the molecule is in one of 24 singly 
charged or 60 doubly charged states with equal probability, 
and decreases to (4 x + 24 x 1 + 60 x 2)/(4 + 24 + 60) = 
18/11 « 1.64 at the second step, when 4 additional neutral 
states become available, etc. Furthermore, we find that each 
Coulomb-blockade step shows a characteristic fine structure, 
which we discuss below. 

The calculation of the differential conductance dl/dV as a 
function of source-drain voltage V and gate voltage V g shows 



the usual Coulomb diamonds; see Fig. [2] Close to the degen- 
eracy points between different charge states we observe a rela- 
tively complex fine structure, corresponding to the steps in the 
inset of Fig.^ We assume very low temperatures, ksT -C J, 
to exhibit the structure more clearly. At higher temperatures 
the peaks in dl/dV are thermally broadened. In the follow- 
ing, we briefly explain the physics behind the fine structure, 
taking Fig.|2ja) as an example. 

Since the Ceo spin Sc eo , the spin of the nitrogen atom, 5n, 
and the total spin S (where S — \Sc eo — Su\, ■■■iSc 60 + $n) 
are good quantum numbers, the exchange energy is 

£cxc = -^[S(S+1)-Sc 60 (Sc 60 + 1)-Sn(Sn + 1)], (13) 

which leads to the level splitting illustrated in Fig. [3] The 
structure in Fig. |2ja) arises from transitions between charge 
states —1 and —2, taking spin excitations into account. In 
equilibrium (V = V) only the ground state of the dot is 
occupied, which is the Ceo state with Sc 60 = 1/2 and 
S = 2 for V g smaller than the degeneracy point Vj 1 and the 
C 60 2 ~ state with S Ceo = 1 and S = 5/2 for V e > V°; 
cf. Fig.|3ja). For V g < V® the first peak in dl/dV at nonzero 
V originates from the transition with Sc 6a = 1/2 — ► 1 
and 5 = 2^5/2, corresponding to a gain of exchange 
energy of AE CXC = — 0.75meV. The second peak results 
from a transition with Sc B0 — 1/2 — > 0, S — 2 — > 3/2, 
and AE CXC = +0.75 meV. Simultaneously, the transitions 
with Sc 60 = 1/2 -> 0, S = 1 -> 3/2, and A£ oxc = 
-1.25 meV and Sc 60 = 1/2 ->■ 1, S = 1 -> 3/2, and 
Ai5 exc = — 0.25meV are enabled [dashed lines in Fig.|3ja)]. 
Although energetically possible, these transitions are not ex- 
cited at lower source-drain voltages, because the lower levels 
are unoccupied. The last two peaks belong to transitions with 
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(a) 



V g (mV) 

FIG. 2: Gray-scale plots of the differential conductance dl/dV as 
a function of source-drain voltage V and gate voltage V g for T — 
O.IK. Shown are two particular ranges of gate voltages close to 
the degeneracy points between charge states —1 and —2 (a),(c) and 
between —2 and —3 (b),(d). (a) and (b) show results for vanishing 
magnetic field and (c) and (d) for B = 2 T. 
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FIG. 3: Energy levels and all allowed transitions between many- 
particle states with one (TV = 1) and two (TV = 2) electrons, taking 
into account spin excitations, (a) Situation with the TV = 1 multiplet 
lower in energy than the TV = 2 multiplet. (b) Reverse case. 



S Ceo = 1/2 -> 1, S = 1 -> 1/2, and AE CXC = +1.25 meV 
and S = 1/2 ->• 1, 5 = 2 -> 3/2, and A£ cxc = +1.75 mcV. 
Note that the values of AE exc account for the level spacing. 

The situation is different for V g significantly larger than V®, 
where we observe only two peaks; cf. Fig. Eta). As soon 
as the transition from the Cqq 2 ^ ground state into the low- 
est C 6 cT state with S Ceo = 1 ->• 1/2, S = 5/2 -» 2, 
and AE CXC = +0.75 meV becomes possible, the transi- 
tions corresponding to AE CXC = — 0.75 me V, — 1.75meV, 
+0.25 me V, and — 1.25meV [dashed lines in Fig. EJb)] are 
also enabled. Again the latter four would be energetically pos- 
sible at lower V, but do not appear as peaks of dl / dV, since 
the corresponding lower levels are unoccupied. In the vicinity 
of V® we find that the slope of several lines abruptly changes 
sign. This corresponds to the situation where two levels con- 
nected in Fig.[3]by a transition cross as V g is varied. The fine 
structure in Fig.|2jb) can be discussed analogously. The struc- 
ture is different for all degeneracy points and can thus serve 
as a fingerprint of the particular charge transition. This should 
be useful since the zero of the V s axis is often shifted signifi- 
cantly from one experiment to the next. 

Selection rules for single-electron tunneling require that 
AS C60 = ±1/2 and AS = ±1/2. The different brightness of 
the peaks in Fig.|2]is correlated with the number of transitions 
that are possible at a given source-drain voltage. Each allowed 
transition may be thought of as one current channel. 

Experimentally, the magnetic origin of the fine structure is 



most conclusively tested by observing the behavior in a mag- 
netic field. For ionized Cgo in lattices and in solution, the 
orbital moment is quenchedi^iS We assume that the fields 
generated by the electrodes in a break junction are also suf- 
ficiently strong to quench the orbital moment. Then the 
molecule couples to a magnetic induction B only through the 
spin moments, described by the new Hamiltonian 

H' = H- gp B BS^ 60 ~ gf i B BS^ =H- gi i B BS z . (14) 

Here, /is is the Bohr magneton and g is the g factor, which 
is g w 2 for both the nitrogen spin Sn and the Ceo spin. We 
choose a many-particle basis of simultaneous eigenstates of 
n-d, Sc 60 , S, and S z . Then the only difference is that addi- 
tional Zeeman energies appear in our expression for the tran- 
sition rates, Eq. dlOl i. In Figs.|2jc) and|2jd) we show dl/dV 
for the same parameters as in Figs. |2ja) and |2jb) but with 
B = 2 T. As expected, the peaks split, but in addition sev- 
eral peaks are absent since they are not allowed by the selec- 
tion rules. For example, for V g > Vj 1 the first peak is due 
to a transition with S Cbo = 1 — ► 1/2, S = 5/2 — > 2, and 
TV = 2 — ► 1. Since the initial state has all spins aligned in 
parallel, one electron tunneling out of the dot can only reduce 
S z so that there is only a single peak in dl/dV. 

To summarize, we have presented a theory for transport 
through a single N@Cgo molecule weakly coupled to metallic 
electrodes. Our results for the differential conductance dl/dV 
as a function of the source-drain and gate voltages show 
Coulomb-blockade and exhibit a characteristic fine structure 
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of the Coulomb-blockade peaks due to the coupling of the Cgo 
spin to the spin of the encapsulated nitrogen atom. 
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